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OBJECT  OF  TASK 

To  improve  existing  knowledge  of  gamma  and  neutron  shielding  properties  of 
shelters. 


ABSTRACT 

The  use  of  the  albedo  approach  for  the  analytical  prediction  of  attenuation 
characteristics  of  two-legged  rectangular  ducts  transmitting  gamma  rays  is  extended 
to  cases  involving  the  asymmetrical  location  of  source  and/or  detector.  Illustrative 
problems  are  presented  and  solved. 


II 


I.  INTRODUCTION 

The  problem  of  radiation  streaming  down  ducts  and  passageways  in  a  shield 
has  many  facets.  The  part  of  the  problem  relating  to  attenuation  of  gamma  radia¬ 
tion  through  two-legged  rectangular  ducts,  with  symmetrical  placement  of  source 
and  detector,  has  been  the  subject  of  several  recent  investigations,  both  experimen¬ 
tal  and  analytical  (References  1-3). 

This  note  carries  on  the  program  of  analysis  initiated  in  Reference  1  and 
undertakes  an  analytical  solution  of  similar  problems  for  cases  where  the  source 
and/or  the  detector  are  not  located  on  the  axis  of  symmetry  of  the  legs,  but  are  to 
one  side. 


II.  BACKGROUND 

Reference  1  has  displayed  in  some  detail  a  method  of  predicting  the 
attenuation  of  gamma  radiation  as  it  passes  down  a  two-legged  rectangular  duct 
within  a  thick  shield.  There  is  therefore  no  need  here  for  a  precise  exposition  of 
the  general  problem  and  the  precise  assumptions  that  are  necessary  to  obtain  ana¬ 
lytical  solutions.  It  is  assumed  that  the  reader  has  become  acquainted  with  such 
matters  from  Reference  1 . 

An  outline  of  the  problem,  however,  is  desirable.  Its  main  aspects  are  as 
follows: 

(a)  We  are  concerned  with  ducts  of  rectangular  cross  section,  having  more 
than  one  leg,  with  the  legs  intersecting  perpendicularly.  These  ducts  may  be  as 
small  as  an  air  shaft,  or  they  may  be  as  large  as  an  entranceway  to  a  personnel 
shelter  designed  to  provide  protection  against  atomic  bomb  radiation. 

(b)  The  minimum  dimension  of  the  ducts  is  several  inches,  and  the  lengths 
are  less  than,  say,  a  hundred  feet.  The  ratio  of  length  to  width  or  height  is  con¬ 
sidered  to  be  three  or  more.  (Reference  1,  however,  indicates  surprisingly  accurate 
results  at  ratios  almost  down  to  unity.) 


(c)  A  high  order  of  accuracy  is  sought  by  eliminating  any  known  approxi¬ 
mations  giving  errors  greater  than  10  percent;  however,  an  overall  ratio  between 
analysis  and  experiment  of  not  over  1 .5  is  considered  adequate. 

(d)  The  analytical  approach  employs  the  "albedo"  concept.  Multiple 
scattering  off  walls  is  ignored,  for  lack  of  adequate  means  of  handling  this  aspect 
of  the  problem  at  present. 

(e)  The  wall  material  is  uniform  in  composition  and  density. 

III.  THE  PROBLEM  OF  OFF-CENTERED  SOURCE  AND/OR  DETECTOR 

In  Reference  I  the  general  problem  was  simplified  by  assuming  all  the  source 
radiation  to  be  concentrated  at  the  center  of  the  enhance  to  the  first  duct  leg  (or 
assuming  by  symmetry  that  it  could  be  so  regarded),  and  by  having  the  detector  at 
the  center  of  the  exit  of  the  second  duct  leg.  There  are  circumstances,  however, 
where  this  assumption  should  not  or  cannot  be  made.  As  a  further  refinement, 
therefore,  one  should  investigate  the  cases  in  which  source  and/or  detector  are 
placed  at  the  "farther"  or  "nearer"  wall,  of  the  duct  system.  Figures  I,  2,  and  3 
illustrate  the  various  situations  studied. 

It  can  be  readily  demonstrated  that  the  position  of  the  source  away  from  the 
center,  in  the  direction  perpendicular  to  the  plane  established  by  the  axes  of  the 
duct  legs,  does  not  significantly  affect  the  analytical  results;  whereas  change  of 
location  in  the  plane  of  the  axis  does.  We  therefore  need  not  concern  ourselves 
about  location  at  the  other  two  sides  of  the  access  openings  (the  floor  and  ceiling 
of  a  horizontal  entranceway,  for  example).  The  "far"  or  "near"  wall  cases  are  to 
be  considered  only. 

The  procedure  in  each  case  is  almost  precisely  the  same  as  that  given  in 
Reference  1,  and  we  shall  not  include  all  detailed  derivations  in  this  treatise.  Only 
the  pertinent  results  will  be  noted. 

Table  I  indicates  how  the  contributions  to  the  detector  response  in  various 
cases  are  readily  established  by  comparison  with  that  of  the  fundamental  case  ana¬ 
lyzed  in  Reference  1,  the  so-called  C-C  case.  By  using  this  table,  one  can  obtain 
and  add  the  various  wall  scattering  contributions  to  the  overall  detector  readings, 
for  the  various  cases.  We  can,  then,  for  each  case  use  the  same  fundamental 
formulas  provided  in  Reference  1  for  the  C-C  case: 
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Ft  =  F,  f2 


F2  =  (4p |  p2  p3)  Gtot 
GTOT  =  Gb  +  Gt1  +  Gt2  +  Gs 

in  which  Fj  is  the  overall  attenuation  factor  D/DQ,  F^  is  the  attenuation  factor  for 
the  first  leg  (not  particularly  sensitive  to  location  of  source  across  the  entranceway)* 
F2  is  the  attenuation  factor  for  the  second  leg  (of  prime  concern  to  us  here),  and  the 
G's  are  terms  or  combinations  of  terms  (largely  geometrical  in  nature)  which  permit 
ready  calculation  of  the  attenuation  factors. 

The  various  G-terms  for  the  several  cases  are  listed  below.  For  G^  the 
following  are  readily  derived. 


Case  C-C: 


1  2  3  a4 

G  =  +  r — .:■* ■!■■■■:  +  - — +  — 4 


1  +  p,  p2(i  +  p2)  1  -  p,  1  -  p2 


Case  F-C: 


G.  = 


*1  u2  u3  4 

+  — —  +  — - —  +  - — 


b  1  +  p,  p2  i-p 


'  o-p2)‘ 


Case  N-C: 


al  a2  °3 

°b  =  ^ +  p2(i  +  p/  +  ^  +  °4 


Case  C-F: 


r  ,  °2  ,  °3  a4 

G.  =  a,  +  r — r - —  +  - -  + 

b  I 


P2  0  +  P2)  (1  -  M 
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Case  F-F: 


«...  ♦  -X *  —  3—  4.  — a_ 

b  (1-P,)2  (1-P2)2 


Case  N-F: 


a2  °3 

G i  -  a,  +  .  1  — „  +  — — — — ■  ■+■  q  : 

b  1  P20  +  P2)2  0  -  P,)2  4 


v 


Case  C-N: 


Gl  = 


ql  q2  q4 

+  p  P2  ^  +  P2^  °3 


Case  F-N: 


°1  °2  a4 

Gu  =  2  +  p2  +  a3  +  .  2 


0+P,) 


0  ~P2) 


Case  N-N: 


1  2 

G  =  - -  +  - -  +  a«  +  a  A 

b  o  +  p,)2  P2o+P2r  3  4 


Gjj  is  the  same  for  cases  C-C,  F-C,  and  N-C;  it  is  as  given  in  Reference  1. 
For  cases  C-F/  F-F/  and  N-F/  it  is  twice  as  much  as  for  the  fundamental  C-C  case. 
For  cases  C-N/  F-N,  and  N-N,  it  is  zero. 

Gj2  >s  the  same  for  cases  C-C,  C-F,  and  C-N;  it  is  as  given  in  Reference  1 . 
For  cases  F-C,  F-F,  and  F-N,  it  is  twice  as  much  as  for  the  fundamental  C-C  case. 
For  cases  N-C,  N-F,  and  N-N,  it  is  zero. 
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Table  I.  Relative  Values  of  Scattering  Areas  and  Associated  Dose  Contributions 


^'s^ase 

Areas\ 

C-C 

F-C 

N-C 

C-F 

F-F 

N-F 

C-N 

F-N 

N-N 

1 

i 

1 

1 

(a) 

(a) 

(a) 

(c) 

(c) 

(c) 

2 

1 

(b) 

(d) 

1 

(b) 

(d) 

1 

(b) 

(d) 

3 

1 

1 

1 

(a) 

(a) 

(a) 

(c) 

(c) 

(c) 

4 

1 

(b) 

(d) 

1 

(b) 

(d) 

1 

(b) 

(d) 

5 

1 

1 

1 

2 

2 

2 

0 

0 

0 

6 

1 

2 

0 

1 

2 

0 

1 

2 

0 

7 

1 

1 

1 

2 

2 

2 

0 

0 

0 

8 

1 

2 

0 

1 

2 

0 

1 

2 

0 

Key  to  cases:  C  indicates  location  in  Center. 

F  indicates  location  at  Farther  Wall. 

N  indicates  location  at  Nearer  Wall. 

The  first  letter  denotes  location  of  source; 
the  second  letter  denotes  location  of  detector. 
For  example,  C-F  means  that  the  source  is 
located  at  the  center  of  the  entrance  to  the 
first  leg,  and  the  detector  is  located  on  the 
farther  wall  at  the  end  of  the  second  leg. 


Key  to  tabulated  entries:  Values  are  relative  to  the  C-C  case,  which  is  given 

in  Tables  I  and  III  of  Reference  1. 


W. 


(a) 

(b) 

(c) 

(d) 


indicates 

indicates 

indicates 

indicates 

W„ 


Notei  pi  =  Ti^  '  p2  = 


2Li 


1  +  P  ]  /  or  1/(1  ~  P |)» 
1  +  P2  f  or  1/(1  -  P2)* 
1  -  Pi  ,  or  1/(1  +P}). 

1  ~  P2  '  or  1/(1  + 
p,  — H-  . 

p3  2  L2 
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The  various  expressions  for  Gs  are  readily  obtained  by  carrying  through/  for 
each  case,  the  same  sort  of  derivation  given  In  Reference  1  for  the  C-C  case.  In 
general,  the  following  expressions  are  obtained. 


,  “1 
tan 

K,W, 

2  L,  (1  -  p2) 

tan  * 

k2w2 

2l20  -  Pi) 

90°  ■ 

•  a\  -  a2 

K3  Z  N  K(0S,  Ec) 

2  L2(l  -  P,)3  (1  -  P2)? 

The  values  of  Kj,  K2,  and  Kg  are  given  in  Table  II  for  the  various  cases. 

Table  II.  The  Values  of  the  Constants  K  in  Equations 


Case 

Ki 

k2 

K3 

C-C 

1 

1 

1 

F-C 

2 

1 

2 

N-C 

- 

- 

0 

C-F 

1 

2 

2 

F-F 

2 

2 

4 

N-F 

- 

- 

0 

C-N 

- 

- 

0 

F-N 

- 

- 

0 

N-N 


0 


IV.  EXAMPLE  NO.  1 


Problem.  Consider  a  two-legged  duct,  of  square  cross  section,  28.2  centi¬ 
meters  on  a  side.  The  walls  are  of  lightweight  concrete,  124  pcf.  Each  leg  is 
80  centimeters  long.  A  point  source,  with  a  photon  energy  assumed  to  be  0.34 
Mev,  is  located  at  the  "far"  side  of  the  entrance  to  the  duct.  The  detector  is 
located  on  the  center  line  at  the  exit.  What  is  the  overall  attenuation  of  the  duct, 
using  the  response  of  a  detector  one  centimeter  away  from  the  source  in  free  space 
as  a  reference? 

Solution.  The  solution  is  carried  out  as  indicated  in  Table  III  herein,  which 
has  essentially  the  same  format  as  prescribed  for  this  kind  of  problem  by  Table  V  in 
Reference  1 .  The  case  herein  is  the  F-C  case  and  the  formulas  must  be  adjusted  as 
described  in  Section  III  above.  The  attenuation  is  1.72  x  10“^. 


V.  EXAMPLE  NO.  2 

Problem.  Figure  4  illustrates  a  duct  system  with  a  parallel,  broad  beam  of 
0.5-Mev  gamma  photons  incident  at  an  angle  ij)  =  38°  with  the  normal  to  the  far 
side  of  the  first  duct  leg.  The  duct  dimensions  and  materials  are  the  same  as  in  the 
previous  example,  except  that  the  length  of  the  first  leg  is  91  centimeters.  If  the 
outside  dose,  within  the  beam,  is  100,000  roentgens,  what  dose  is  registered  by  the 
detector  at  the  end  of  the  second  leg?  (Assume  no  direct  penetration  of  the  shield 
itself.) 

Discussion.  Let  us  briefly  review  the  albedo  concept  before  proceeding  with 
the  solution.  (Figure  5,  which  is  the  same  as  Figure  1  of  Reference  1,  Is  repro¬ 
duced  for  reference  purposes.) 

The  basic  albedo  equation.  Equation  1  of  Reference  1,  can  be  written  as 
follows: 


D  =  (D0)  (-L)  (A  a  cos9,)  (-L-) 

ri  r2 
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Let  u$  consider  the  factors  on  the  right-hand  side  of  the  equatton,  from  left  to 
right.  The  first  term  gives  the  reading  of  the  detector  at  the  original  reference 
point.  Multiplied  by  the  next  factor,  it  provides  a  detector  response  at  the  scat¬ 
tering  area  without  scattering,  and  thus  indicates  the  level  of  radiation  striking 
the  surface.  Multiplied  then  by  the  third  factor,  it  gives  the  detector  reading  a 
unit  distance  away  from  the  scattering  surface  after  and  Including  only  the  scat¬ 
tering  from  the  area  A.  (This  latter  detector  reading  is  somewhat  hypothetical 
unle  the  dimensions  of  the  scattering  surface  are  much  less  than  unity  —  it  is 
nevertheless  a  useful  "mental  crutch.")  When  the  final  factor  is  also  included,  we 
get  the  detector  response  at  the  distance  ^  from  the  scattering  surface.  Thus  the 
scattering  surface,  if  its  dimensions  are  small  compered  to  ^  acts  essentially  as  a 
point  source  for  radiation  going  along  path  r2* 

Under  conditions  where  the  incident  radiation  does  not  proceed  from  a  point 
source  but  is  a  broad,  parallel  beam  incident  at  angle  ip  to  the  normal,  it  is  readily 
seen  that  the  first  two  factors  are  replaced  simply  by  the  detector  reading  in  the 
unscattered  beam,  which  we  will  call  D|J.  Thus  the  albedo  equation  would  be: 

D  =  A  a  cos  ip  f - 

V 

In  such  case,  we  see,  the  scattering  area  A  acts  as  If  It  were  a  point  source  with 
sufficient  strength  for  the  detector  reading  one  unit  distance  from  it  to  give  a 
value  of  D|J  A  a  cos  ip. 

Solution.  The  problem  quickly  reduces  to  the  sort  of  two-legged  duct  problem 
we  have  discussed  in  this  paper,  the  F-C  case,  if  we  take  the  mid-point  of  the 
scattering  area  at  the  entrance  side  wall  as  the  location  of  the  source  in  the 
reduced  problem.  The  strength  of  the  source  in  the  reduced  problem,  which  we  will 
call  Dq,  is,  as  explained  in  the  above  discussion: 


I  ^ 

Dq  =  Dy  Aa  cos  ip 

=  Dy  (Wj  H  tani/i)  a  cos  ip 
=  Dy  W|  H  a  sin^ 
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Table  III.  Solution  to  Problem  No.  1/  continued  (page  two  of  two) 


For  the  F-C  case: 

a1  a2  .  a3  °4 

b  1  +  P 1  ^2  1  "h  (1  -  P2)2 

=  0.039  +  0.085  +  0.056  +  0.068 

Gb  =  0.248 

(1  "Pi)  a5  "  2  a7 

Gt  ~  ,  9  3 

2m;l2o-p,)2  o-p2)3 

(1  -  p2)  a6  -  2  p2ag 

f  2  2 
4aL1P2<1“P|)  0“P2) 

0.1214  0.0286 

3.60  0.386 

=  0.0337  +  0.0740 

Gt  =  0.108 

’3  3) 

\(>  -Pi)3  0-P2)3 

=  3.78/25.0 

Gs  =  0.151 

GTOT  =  0,507 

F1  =  1/L12 

F2  "  4Pl  ^2  ^3  gTOT 

Ft  =  f,  f2 

F1  =  1.10  x  10"7 
F2  =  1.56  x  10-4 
Ft  =  1.72  x  10"6 
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Figure  4.  Illustration  for  Example  No. 
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Figure  5.  Scattering 


Substituting  the  ptven  numerical  data  and  obtaining  the  value  of  a  from  Figure  4, 
Reference  1,  we  get: 


Dq  =  (100,000)  (28.2)  (28.2)  (0.616)  (0.0148) 

=  725, 000  roentgens-cm  . 

(The  fact  that  this  figure  is,  astonishingly,  larger  than  the  value  of  Dy,  the  reading 
in  open  air.  Is  due  to  the  hypothetical  nature  of  the  concept  that  the  radiation 
scattered  from  an  area  of  the  side  wall  near  the  opening  is  considered  to  be  concen¬ 
trated  into  a  single  point.) 

The  radiation  scattered  off  the  far  wall  at  the  entrance  may  be  considered, 
without  undue  error,  as  coming  from  a  single  point  on  the  wall,  at  a  distance  from 
the  opening  equal  to  1/2  Wj  tan^>,  which  here  has  a  numerical  value  of  11  centi¬ 
meters.  For  the  reduced  problem,  then,  the  length  of  the  first  leg  is  91  minus  11, 
or  80  centimeters. 

The  next  quantity  needed  is  the  representative  energy  for  the  gamma  photons 
after  scattering.  Assuming  that  most  of  the  albedo  is  due  to  single  scattering,  we 
can  use  the  usual  Compton  scattering  relationship  (Reference  5,  p.  146)  to  deter¬ 
mine  the  scattered  photon  energy.  From  geometrical  relationships,  the  average 
scattering  angle  is  found  to  be  62  degrees;  by  use  of  the  Compton  relationship,  the 
scattered  photon  energy  is  found  to  be  about  0.34  Mev. 

At  this  point  we  solve  the  reduced  problem  for  a  two-legged  duct  attenuation, 
using  the  F-C  case.  For  this  particular  example,  the  solution  has  already  been  , 
obtained  in  Example  1,  above.  The  dose  recorded  by  the  detector  at  the  end  of  the 
second  leg  of  the  duct  is  therefore  found  to  be  (1 .72  x  10"^)  x  (0.725  x  10  ),  or 
1 .25  roentgens. 


VI.  CONCLUDING  REMARKS 

We  have  no  experimental  data  to  confirm  adequately  the  type  of  analysis 
derived  and  utilized  herein.  The  results  of  this  work  therefore  cannot  be  said  to 
have  achieved  firm  scientific  status.  It  Is,  however,  based  on  a  plausible  exten 
slon  of  principles  and  techniques  which  are  more  firmly  established,  and  for  the 
time  being  will  serve  to  permit  engineering  calculations  for  these  and  related 
practical  situations.  The  need  for  confirming  experimental  work  is  clear. 
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